Extended Cleavage Specificity of two Hematopoietic Serine Proteases from a Ray-Finned Fish, the Spotted Gar (Lepisosteus oculatus).

The extended cleavage specificities of two hematopoietic serine proteases originating from the ray-finned fish, the spotted gar (Lepisosteus oculatus), have been characterized using substrate phage display. The preference for particular amino acids at and surrounding the cleavage site was further validated using a panel of recombinant substrates. For one of the enzymes, the gar granzyme G, a strict preference for the aromatic amino acid Tyr was observed at the cleavable P1 position. Using a set of recombinant substrates showed that the gar granzyme G had a high selectivity for Tyr but a lower activity for cleaving after Phe but not after Trp. Instead, the second enzyme, gar DDN1, showed a high preference for Leu in the P1 position of substrates. This latter enzyme also showed a high preference for Pro in the P2 position and Arg in both P4 and P5 positions. The selectivity for the two Arg residues in positions P4 and P5 suggests a highly specific substrate selectivity of this enzyme. The screening of the gar proteome with the consensus sequences obtained by substrate phage display for these two proteases resulted in a very diverse set of potential targets. Due to this diversity, a clear candidate for a specific immune function of these two enzymes cannot yet be identified. Antisera developed against the recombinant gar enzymes were used to study their tissue distribution. Tissue sections from juvenile fish showed the expression of both proteases in cells in Peyer's patch-like structures in the intestinal region, indicating they may be expressed in T or NK cells. However, due to the lack of antibodies to specific surface markers in the gar, it has not been possible to specify the exact cellular origin. A marked difference in abundance was observed for the two proteases where gar DDN1 was expressed at higher levels than gar granzyme G. However, both appear to be expressed in the same or similar cells, having a lymphocyte-like appearance.

Functional Divergence in Solute Permeability between Ray-Finned Fish-Specific Paralogs of aqp10.

Aquaporin (Aqp) 10 is a member of the aquaglyceroporin subfamily of water channels, and human Aqp10 is permeable to solutes such as glycerol, urea, and boric acid. Tetrapods have a single aqp10 gene, whereas ray-finned fishes have paralogs of this gene through tandem duplication, whole-genome duplication, and subsequent deletion. A previous study on Aqps in the Japanese pufferfish Takifugu rubripes showed that one pufferfish paralog, Aqp10.2b, was permeable to water and glycerol, but not to urea and boric acid. To understand the functional differences of Aqp10s between humans and pufferfish from an evolutionary perspective, we analyzed Aqp10s from an amphibian (Xenopus laevis) and a lobe-finned fish (Protopterus annectens) and Aqp10.1 and Aqp10.2 from several ray-finned fishes (Polypterus senegalus, Lepisosteus oculatus, Danio rerio, and Clupea pallasii). The expression of tetrapod and lobe-finned fish Aqp10s and Aqp10.1-derived Aqps in ray-finned fishes in Xenopus oocytes increased the membrane permeabilities to water, glycerol, urea, and boric acid. In contrast, Aqp10.2-derived Aqps in ray-finned fishes increased water and glycerol permeabilities, whereas those of urea and boric acid were much weaker than those of Aqp10.1-derived Aqps. These results indicate that water, glycerol, urea, and boric acid permeabilities are plesiomorphic activities of Aqp10s and that the ray-finned fish-specific Aqp10.2 paralogs have secondarily reduced or lost urea and boric acid permeability.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes.

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Conserved enhancers control notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Suppressors of cGAS-STING are downregulated during fin-limb regeneration and aging in aquatic vertebrates.

During the early stages of limb and fin regeneration in aquatic vertebrates (i.e., fishes and amphibians), blastema undergo transcriptional rewiring of innate immune signaling pathways to promote immune cell recruitment. In mammals, a fundamental component of innate immune signaling is the cytosolic DNA sensing pathway, cGAS-STING. However, to what extent the cGAS-STING pathway influences regeneration in aquatic anamniotes is unknown. In jawed vertebrates, negative regulation of cGAS-STING activity is accomplished by suppressors of cytosolic DNA such as Trex1, Pml, and PML-like exon 9 (Plex9) exonucleases. Here, we examine the expression of these suppressors of cGAS-STING, as well as inflammatory genes and cGAS activity during caudal fin and limb regeneration using the spotted gar (Lepisosteus oculatus) and axolotl (Ambystoma mexicanum) model species, and during age-related senescence in zebrafish (Danio rerio). In the regenerative blastema of wounded gar and axolotl, we observe increased inflammatory gene expression, including interferon genes and interleukins 6 and 8. We also observed a decrease in axolotl Trex1 and gar pml expression during the early phases of wound healing which correlates with a dramatic increase in cGAS activity. In contrast, the plex9.1 gene does not change in expression during wound healing in gar. However, we observed decreased expression of plex9.1 in the senescing cardiac tissue of aged zebrafish, where 2'3'-cGAMP levels are elevated. Finally, we demonstrate a similar pattern of Trex1, pml, and plex9.1 gene regulation across species in response to exogenous 2'3'-cGAMP. Thus, during the early stages of limb-fin regeneration, Pml, Trex1, and Plex9.1 exonucleases are downregulated, presumably to allow an evolutionarily ancient cGAS-STING activity to promote inflammation and the recruitment of immune cells.

Ancient vertebrate dermal armor evolved from trunk neural crest.

Bone is an evolutionary novelty of vertebrates, likely to have first emerged as part of ancestral dermal armor that consisted of osteogenic and odontogenic components. Whether these early vertebrate structures arose from mesoderm or neural crest cells has been a matter of considerable debate. To examine the developmental origin of the bony part of the dermal armor, we have performed in vivo lineage tracing in the sterlet sturgeon, a representative of nonteleost ray-finned fish that has retained an extensive postcranial dermal skeleton. The results definitively show that sterlet trunk neural crest cells give rise to osteoblasts of the scutes. Transcriptional profiling further reveals neural crest gene signature in sterlet scutes as well as bichir scales. Finally, histological and microCT analyses of ray-finned fish dermal armor show that their scales and scutes are formed by bone, dentin, and hypermineralized covering tissues, in various combinations, that resemble those of the first armored vertebrates. Taken together, our results support a primitive skeletogenic role for the neural crest along the entire body axis, that was later progressively restricted to the cranial region during vertebrate evolution. Thus, the neural crest was a crucial evolutionary innovation driving the origin and diversification of dermal armor along the entire body axis.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

PML and PML-like exonucleases restrict retrotransposons in jawed vertebrates.

We have uncovered a role for the promyelocytic leukemia (PML) gene and novel PML-like DEDDh exonucleases in the maintenance of genome stability through the restriction of LINE-1 (L1) retrotransposition in jawed vertebrates. Although the mammalian PML protein forms nuclear bodies, we found that the spotted gar PML ortholog and related proteins in fish function as cytoplasmic DEDDh exonucleases. In contrast, PML proteins from amniote species localized both to the cytoplasm and formed nuclear bodies. We also identified the PML-like exon 9 (Plex9) genes in teleost fishes that encode exonucleases. Plex9 proteins resemble TREX1 but are unique from the TREX family and share homology to gar PML. We also characterized the molecular evolution of TREX1 and the first non-mammalian TREX1 homologs in axolotl. In an example of convergent evolution and akin to TREX1, gar PML and zebrafish Plex9 proteins suppressed L1 retrotransposition and could complement TREX1 knockout in mammalian cells. Following export to the cytoplasm, the human PML-I isoform also restricted L1 through its conserved C-terminus by enhancing ORF1p degradation through the ubiquitin-proteasome system. Thus, PML first emerged as a cytoplasmic suppressor of retroelements, and this function is retained in amniotes despite its new role in the assembly of nuclear bodies.

Genome structures resolve the early diversification of teleost fishes.

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

Na+/Cl- cotransporter 2 is not fish-specific and is widely found in amphibians, non-avian reptiles, and select mammals.

Solute carrier 12 (Slc12) is a family of electroneutral cation-coupled chloride (Cl-) cotransporters. Na+/K+/2Cl- (Nkcc) and Na+/Cl- cotransporters (Ncc) belong to the Nkcc/Ncc subfamily. Human and mouse possess one gene for the Na+/Cl- cotransporter (ncc gene: slc12a3), whereas teleost fishes possess multiple ncc genes, slc12a3 (ncc1) and slc12a10 (ncc2), in addition to their species-specific paralogs. Amphibians and squamates have two ncc genes: slc12a3 (ncc1) and ncc3. However, the evolutionary relationship between slc12a10 and ncc3 remains unresolved, and the presence of slc12a10 (ncc2) in mammals has not been clarified. Synteny and phylogenetic analyses of vertebrate genome databases showed that ncc3 is the ortholog of slc12a10, and slc12a10 is present in most ray-finned fishes, coelacanths, amphibians, reptiles, and a few mammals (e.g., platypus and horse) but pseudogenized or deleted in birds, most mammals, and some ray-finned fishes (pufferfishes). This shows that slc12a10 is widely present among bony vertebrates and pseudogenized or deleted independently in multiple lineages. Notably, as compared with some fish that show varied slc12a10 tissue expression profile, spotted gar, African clawed frog, red-eared slider turtle, and horse express slc12a10 in the ovaries or premature gonads. In horse tissues, an unexpectedly large number of splicing variants for Slc12a10 have been cloned, many of which encode truncated forms of Slc12a10, suggesting that the functional constraints of horse slc12a10 are weakened, which may be in the process of becoming a pseudogene. Our results elaborate on the evolution of Nkcc/Ncc subfamily of Slc12 in vertebrates.NEW & NOTEWORTHY slc12a10 is not a fish-specific gene and is present in a few mammals (e.g., platypus and horse), non-avian reptiles, amphibians, but was pseudogenized or deleted in most mammals (e.g., human, mouse, cat, cow, and rhinoceros), birds, and some ray-finned fishes (pufferfishes).

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development.

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Genome of the Rio Pearlfish (Nematolebias whitei), a bi-annual killifish model for Eco-Evo-Devo in extreme environments.

The Rio Pearlfish, Nematolebias whitei, is a bi-annual killifish species inhabiting seasonal pools in the Rio de Janeiro region of Brazil that dry twice per year. Embryos enter dormant diapause stages in the soil, waiting for the inundation of the habitat which triggers hatching and commencement of a new life cycle. Rio Pearlfish represents a convergent, independent origin of annualism from other emerging killifish model species. While some transcriptomic datasets are available for Rio Pearlfish, thus far, a sequenced genome has been unavailable. Here, we present a high quality, 1.2 Gb chromosome-level genome assembly, genome annotations, and a comparative genomic investigation of the Rio Pearlfish as representative of a vertebrate clade that evolved environmentally cued hatching. We show conservation of 3D genome structure across teleost fish evolution, developmental stages, tissues, and cell types. Our analysis of mobile DNA shows that Rio Pearlfish, like other annual killifishes, possesses an expanded transposable element profile with implications for rapid aging and adaptation to harsh conditions. We use the Rio Pearlfish genome to identify its hatching enzyme gene repertoire and the location of the hatching gland, a key first step in understanding the developmental genetic control of hatching. The Rio Pearlfish genome expands the comparative genomic toolkit available to study convergent origins of seasonal life histories, diapause, and rapid aging phenotypes. We present the first set of genomic resources for this emerging model organism, critical for future functional genetic, and multiomic explorations of "Eco-Evo-Devo" phenotypes of resilience and adaptation to extreme environments.

Convergent losses of SCPP genes and ganoid scales among non-teleost actinopterygians.

Various secretory calcium-binding phosphoprotein (SCPP) genes are expressed in the skin and jaw during the formation of bone, teeth, and scales in osteichthyans (bony vertebrates). Among these mineralized skeletal units is the ganoid scale, found in many fossil actinopterygians (ray-finned fish) but confirmed only in Polypteriformes (bichirs, reedfish) and Lepisosteiformes (gars) among extant clades. Here, we examined SCPP genes in the genome of seven non-teleost actinopterygian species that possess or do not possess ganoid scales. As a result, 39-43 SCPP genes were identified in Polypteriformes and Lepisosteiformes, whereas 22-24 SCPP genes were found in Acipenseriformes (sturgeons, paddlefish) and Amiiformes (bowfin). Most of these genes form two clusters in the genome of Polypteriformes, Lepisosteiformes, and Amiiformes, and these two clusters are duplicated in Acipenseriformes. Despite their distant phylogenetic relationship, Polypteriformes and Lepisosteiformes retain many orthologous SCPP genes. These results imply that common ancestors of extant actinopterygians possessed a large repertoire of SCPP genes, and that many SCPP genes were lost independently in Acipenseriformes and Amiiformes. Notably, most SCPP genes originally located in one of the two SCPP gene clusters are retained in Polypteriformes and Lepisosteiformes but were secondarily lost in Acipenseriformes and Amiiformes. In Lepisosteiformes, orthologs of these lost genes show high or detectable expression levels in the skin but not in the jaw. We thus hypothesize that many SCPP genes located in this cluster are involved in the formation of ganoid scales in Polypteriformes and Lepisosteiformes, and that their orthologs and ganoid scales were convergently lost in Acipenseriformes and Amiiformes.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system.

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

The bowfin genome illuminates the developmental evolution of ray-finned fishes.

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

Bilateral visual projections exist in non-teleost bony fish and predate the emergence of tetrapods.

In most vertebrates, camera-style eyes contain retinal ganglion cell neurons that project to visual centers on both sides of the brain. However, in fish, ganglion cells were thought to innervate only the contralateral side, suggesting that bilateral visual projections appeared in tetrapods. Here we show that bilateral visual projections exist in non-teleost fishes and that the appearance of ipsilateral projections does not correlate with terrestrial transition or predatory behavior. We also report that the developmental program that specifies visual system laterality differs between fishes and mammals, as the Zic2 transcription factor, which specifies ipsilateral retinal ganglion cells in tetrapods, appears to be absent from fish ganglion cells. However, overexpression of human ZIC2 induces ipsilateral visual projections in zebrafish. Therefore, the existence of bilateral visual projections likely preceded the emergence of binocular vision in tetrapods.

Cellular mechanisms of frontal bone development in spotted gar (Lepisosteus oculatus).

BACKGROUND: The cellular and molecular mechanisms initiating vertebrate cranial dermal bone formation is a conundrum in evolutionary and developmental biology. Decades of studies have determined the developmental processes of cranial dermal bones in various vertebrates and identified possible inducers of dermal bone. However, evolutionarily derived characters of current experimental model organisms, such as non-homologous frontal bones between teleosts and sarcopterygians, hinder investigations of ancestral and conserved mechanisms of vertebrate cranial dermal bone induction. Thus, investigating such mechanisms with animals diverging at evolutionarily informative phylogenetic nodes is imperative. RESULTS: We investigated the cellular foundations of skull frontal bone formation in the spotted gar Lepisosteus oculatus, a basally branching non-teleost actinopterygian. Whole-mount bone and cartilage staining and hematoxylin-eosin section staining revealed that mesenchymal cell condensations in the frontal bone of spotted gar develop in close association with the underlying cartilage. We also identified novel aspects of frontal bone formation: enrichment of F-actin, cellular membranes, and E-cadherin in condensing cells, and extension of podia-like structures from osteoblasts to the frontal bone, which may be responsible for bone mineral transport. CONCLUSION: This study highlights the process of frontal bone formation with dynamic architectural changes of mesenchymal cells in spotted gar, an emerging non-teleost fish model system, illuminating supposedly ancestral and likely conserved developmental mechanisms of skull bone formation among vertebrates.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data.

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.